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ABSTRACT 

We present a high resolution radio image of SDSS 0836+0054 identified recently as the most 
distant radio-detected quasar at a redshift of z = 5.82. The observation was carried out with 
ten antennas of the European VLBI Network, spread from Europe to China and South Africa, 
at 1.6 GHz frequency on 2002 June 8. The source is detected with a total flux density of 
1.1 mJy, equal to its flux density measured in the VLA FIRST survey. We found no indication 
of multiple images produced by gravitational lensing. The radio structure of the quasar at ~ 
10-mas angular resolution appears somewhat resolved. It resembles the radio structure typical 
for lower redshift radio-loud active galactic nuclei. We obtained so far the best astrometric 
position of the source with an accuracy better than 8 mas, limited mainly by the structural 
effects in the phase-reference calibrator source. 
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observations 



1 INTRODUCTION 

The object SDSSp J083643.85+005453.3 (hereafter: SDSS 
0836+0054) is the highest redshift radio-detected quasar known to 
date (z = 5.82) discovered using mul ticolor i maging data from 
\ the Sloan Digital Sky Survey (SDSS) <Fan et aljl200lh . The ob- 
ject is one of the seven z > 5.7 quasars known. Six of them 
form a complete colour-selected flux-limited sample from an area 
of 2870 deg 2 iFan et all20ol . The quasar SDSS 0836+0054 was 
identified with an unresolved radio source within 1" in the VLA 
FIRST survey IWhite et al j 19971) , where the source total flux den- 
sity was 1.11 ±0.15 mJy at 1.4 GHz (21 cm wavelen gth). The total 
flux d ensity in the NRAO VLA Sky Survey (NVSS, ICondon et all 
Il998h 2 was 2.5 ± 0.5 mJy at the same fre quency. 

According to the recent studies (e.g. lBarkanal2002t and ref- 
erences therein), z ~ 6 is close to the reionization epoch of the 
Universe. In the optical spectrum of SDSS 0836+0054, however, 
the Gunn-Peterson effect is not seen. There is detectable emission 
blueward of the Lya line, i.e. the intergalactic medium is already 
highly ionized at z ~ 5.8 along this line of sight (Fan et al. 2001). 

The estimated mass of the central black hole in SDSS 
0836+0054 is 4.8 x 1O 9 M , assuming that the quasar is emit- 
ting at the Eddington luminosity and its apparent flux is not sig- 
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nifica ntly magnified by beaming or gravitational lensing ( F an et alJ 
l200ll) . There are, however, indications that the latter assumptions 
may not be correct. Simulations show that up to one third of z ~ 6 
quasars may be magnified due to gravitational lensing by a fac- 
tor of ten or more. Therefore their black hole masses are overesti- 
mated by the same factor iWvithe & Loebl2002Jbl) . The probabil- 
ity estimations strongly depend on the quasar luminosity function 
used. C onsidering a wider range of possible luminosity functions, 
IComerford et alJ l2002t) conclude that present data allow lensing 
probability of essentially 100%. In turn, the luminosity function 
can be constrained using the detections or absence of lensing events 
with multiple images of z ~ 6 quasars. 

Black hole mass values of z ^ 6 quasars have great impor- 
tance for hierarchical structure formation models in the early Uni- 
verse, placing important constraints on parameters such as the seed 
mass, velocity dispersion, radiative efficiency and accretion lumi- 
nosity. The time needed for a supermassive ~ 10 9 AfQ black hole 
to grow from a stellar-mass seed may be com parable to the age of 
the Universe at this redshift jHaiman & Loebl200ll) . 

The first four z > 5.7 SDSS quasars iFan et all 

1200 ll) including SPSS 0836+0054 are detected in X-rays 
I Brandt et al. 2002). Short exploratory Chandra observations re- 
vealed that SPSS 0836+0054 is consistent with a point-like X- 
ray source iBrandt et all2002tlMathur et all2002llSchwartzl2002t 
Bech told etall2 003). Its position agrees with the optical position 
within the astrometric precision of Chandra (~ 1"). Based on the 
limited data available, the broad-band X-ray and optical properties 
of the extremely distant quasars are apparently similar to those of 
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lower redshift ones , implying no strong evoluti on out to z ~ 6 
iMathur et all2002l) . On the contrary. iBechtold et alJ <2003l) found 
that high-redshift quasars are significantly more X-ray quiet than 
their low-redshift counterparts. 

Any indication of a milli-arcsecond scale "core-jet" radio 
structure revealed by Very Long Baseline Interferometry (VLBI) 
could be considered as a strong support for the existence of a su- 
permassive black hole powering SDSS 0836+0054. These struc- 
tures are common in radio-loud active galactic nuclei (AGN) in 
the broade st range of redshifts from z ~ 0.01 to z ~ 4 
(e.g. IParaei et alj fl99i iFomalont eTall EoOO). According to the 
AGN paradigm, the radio emission originates from incoherent syn- 
chrotron emission in compact (parsec-scale) jets in the close vicin- 
ity of the central black hole. 

Here we describe a 1.6-GHz European VLBI Network (EVN) 
observation of SDSS 0836+0054 and present the high resolution 
radio image of the quasar. Based on our phase-reference VLBI ob- 
servation, we derive the accurate astrometric position of the ob- 
ject. We conclude that the image of the quasar is unlikely to be 
split by gravitational lensing. Throughout this letter, we adopt a 
flat cosmological model with Q m = 0.3, Aa = 0.7 and Ho = 
65 km s _1 Mpc -1 . In this model, 1 milli-arcsecond (mas) angu- 
lar separation corresponds to a linear separation of 6.25 pc at the 
distance of SDSS 0836+0054. 



2 OBSERVATION AND DATA REDUCTION 

We used ten antennas of the EVN to observe SDSS 0836+0054 at 
1.6 GHz (18 cm wavelength) in phase-reference mode. The radio 
telescopes and their parameters are listed in Table Q Due to the 
redshift, the observing frequency corresponds to 1 1 GHz in the rest 
frame of the source. 

Phase-referencing is a technique of extending the signal coher- 
ence time in order to increase the sensitivity of the VLBI array. This 
is done by regularly interleaving observations between the weak 
target source and a nearby strong reference source. Delay, delay 
rate and phase solutions obtained for the phase-reference calibrator 
are interpolated and applied for the target source, thus removing 
most of the phase errors introduced by geometric, ionospheric, at - 
mospheric and instrumental effects (e.g. lBeaslev & Conwavll995T) . 
The ability of the EVN to detect weak (sub-mJy), compact radio 
sources has been dem onstrated recently (e.g. iGarrett et alj|2001t 
iGarrington et alj200ll) . 

The observation of SDSS 0836+0054 took place on 2002 
June 8. We selec ted PKS 0837+012 (J0839+0104; z = 1.123, 
IOwenetalJll995l) as the phase-reference calibrator. The angular 
separation between the target and the reference source is ~ 0.8°. 
The ICRF coordinates of PKS 0837+012 are available from the 
VLBA Calibrator Survey 3 iBeaslev et alJl20o"3> . with an accuracy 
of 0.35 and 0.58 mas in right ascension and declination, respec- 
tively. In principle it allows us to make a precise determination of 
the poorly known position of SDSS 0836+0054 through differential 
astrometry inherently provided by the phase-referencing technique. 
The source J0825+0309 was used as a fringe-finder. 

At the frequency of 1.6 GHz, variations in the ionospheric 
electron column density dominate the propagation effects affect- 
ing the quality of phase transfer between the target and calibra- 
tor sources. These variations have static and dynamic compo- 
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Table 1. EVN telescopes and their characteristics at 1.6 GHz 



Radio telescope (Country) 


Diameter 
(m) 


SEFD a 

(Jy) 


Effelsberg (Germany) 


100 


19 


Hartebeesthoek (South Africa) 


26 


450 


Jodrell Bank Mk2 (United Kingdom) 


25 


320 


Medicina (Italy) 


32 


582 


Nanshan (P.R. China) 


25 


1068 


Noto (Italy) 


32 


784 


Onsala (Sweden) 


25 


390 


Sheshan (P.R. China) 


25 


1130 


Toruh (Poland) 


32 


230 


Westerbork (The Netherlands) 


93 6 


30 



a System Equivalent Flux Density 



equivalent diameter of the phased array 

nents, both spatial and temporal, which limit the maximum achiev- 
ab le dynamic range (D) of VLBI phase-referencing. According 
to iBeaslev & Con way 1 1995), the upper limit of D is given by 
D ~ SOuO" 1 (A/) -1 , where AI is peak daily residual total elec- 
tron content (TEC, in units of 10 17 m~ 2 ), v is the observing fre- 
quency in GHz, 8 is the angular distance between the target and ref- 
erence sources (in degrees). Our observation took place relatively 
close after the peak of the solar cycle. Considering the worst case of 
the maximum TEC as A7 = 5 x 10 ~ 17 m~ 2 (which correspon ds to 
the solar maximum daytime value, IBeaslev & Conwavlll995l) . the 
conservative estimate of the maximum achievable dynamic range 
is D ~ 12. The dynamic component is caused essentially by 
medium-scale travelling ionosphere disturbances (MSTIDs). These 
are time- variable and highly unpredictable. MSTIDs are known to 
vary on the time scale of tens o f minutes to hours and linear scal e 
of hundreds of kilometres (e.g. Ivan Velthoven & Spoelstra|ll992h . 
At the maximum, which normally lasts for up to 5% in time, the 
amplitude of MSTIDs could worsen the upper limit of the achiev- 
able dynamic range by a factor of ~ 2 at our target-reference 
separation comparing to the worst case of static component (e.g. 
IConwav & Vermeulenll993T) . By splitting up the data set into four 
equal time intervals, we did not notice any significant difference 
between our results on this time scale. 

The total observing time was 5 hours with dual circular po- 
larization at all antennas but Sheshan where only left circular po- 
larization was observed. The data were recorded with the Mark IV 
VLBI data acquisition system at the rate of 256 Mbit/s per station 
with 2-bit sampling, resulting in 64 MHz total bandwidth. The cor- 
relation took place at the EVN Data Processor at the Joint Institute 
for VLBI in Europe (JIVE), Dwingeloo, The Netherlands. 

We used 5-min switching cycles between SDSS 0836+0054 
and the phase-reference calibrator quasar J0839+0104. The target 
source was observed for ~ 210-s intervals in each cycle. The total 
on-source time for SDSS 0836+0054 was 2.9 hours. The achievable 
1-er image noise was estimated to be 27 /ijy/beam. 

The NRAO AIPS package (e.g. lDiamonJ 19951) was used for 
initial data calibration and imaging. A priori calibration of visibil- 
ity amplitudes was done using measured system temperatures for 
all antennas. The data for the calibrator and fringe-finder sources 
(J0839+0104 and J0825+0309, respectively) were fringe-fitted in 
AIPS using 3-min solution intervals. The solutions were interpo- 
lated and applied to the data of SDSS 0836+0054 as well. We then 
imaged the phase-reference calibrator source with the conventional 
self-calibration and CLEAN technique (e.g. Walk eJl995t) using the 
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Figure 1. 1.6-GHz VLBI images of the phase-reference calibrator PKS 0837+012 (J0839+0104) with (a) uniform and (b) natural weighting. The positive 
contour levels increase by a factor of %/2- (a) The first contours are drawn at —4.5 and 4.5 mJy/beam. The peak brightness is 112 mJy/beam. The restoring 
beam is 3.5 mas X 2.8 mas at a position angle PA=— 87°. (b) The first contours are drawn at —1.3 and 1.3 mJy/beam. The peak brightness is 234 mJy/beam. 
The restoring beam is 11.6 mas X 9.2 mas at PA=74°. 



AIPS tasks CALIB and IMAGR. The images obtained with two dif- 
ferent data weighting schemes are shown in Fig.Q 

In the uniform weighting scheme, the weights associated with 
the data points are inversely proportional to the local density of vis- 
ibility data, i.e. the longer projected baselines are generally over- 
weighted. As a result, the angular resolution becomes relatively 
high but the image noise level increases. On the other hand, nat- 
ural weighting minimizes the image noise by associating weights 
proportional to the power of —2 of the RMS noise with each data 
point. 

The residual amplitude and phase corrections derived during 
the hybrid mapping of the reference source using natural weighting 
were also applied to the data of the target quasar. The naturally 
weighted image of SDSS 0836+0054 after one cycle of CLEAN 
iterations is shown in Fig. [2] The measured off-source RMS noise 
in the image is 34 ^Jy/beam, close to the theoretical value. 
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3 RESULTS AND DISCUSSION 

The quasar SDSS 0836+0054 is clearly detected with VLBI and ap- 
pears somewhat resolved at a linear resolution of ~ 70 pc (Fig.|2j. 
However, the apparent structural extension may be due to the re- 
maining residual phase errors after phase-reference calibration. 

As seen in Fig. the position of SPSS 0836 +0054 differs 
from the optical position reported bv lFan et a l. ( 2001 ), but is within 
the given error of 200 mas in both right ascension and declination. 
The optical position was used in the correlation of our VLBI data as 
the phase center. Also within the errors, the coordinates are consis- 
tent with the values taken from the VLA FIRST survey. Our values 
OJ2000 = 8 h 36 m 43!8606 and (Wo = 0°54'53'.'232 are deter- 
mined with respect to the accurately known position of the phase- 
reference calibrator J0839+0104. 

However, the differential astrometry is complicated by the fact 
that the structure of our reference source, J0839+0104, is rather 
complex. As seen in Fig.0), there is a low surface brightness ex- 
tended feature to the SW of the brightest component. A similar 
structure in terms of position angle is seen in snapshot VLBA im- 



Figure 2. 1.6-GHz VLBI image of SDSS 0836+0054 with natural weight- 
ing. The positive contour levels increase by a factor of \f 7 l. The first 
contours are drawn at —75 and 75 /^Jy/beam. The peak brightness is 
770 /ijy/beam. The restoring beam is 11.9 mas X 9.4 mas at PA=74°. The 
coordinates are relative to the phase center coincident with the optical posi- 
tion. 



ages taken at higher freque ncies (2.3 and 8.4 GHz) in the VLBA 
Calibrator Survey iBeaslev et all2002T) . With a higher E-W resolu- 
tion (uniform weighting, Fig.^), the "core" splits into two distinct 
components 6.6 mas apart. Supposing that the jet is one-sided, we 
suspect that the true flat-spectrum VLBI core of the quasar is the 
eastern one which is weaker at 1.6 GHz. Most likely this is the 
reference point of which the accurate coordinates are determined 
using 8.4-GHz VLBI observations. On the other hand, the bright- 
est feature in our image may well be a steep-spectrum jet compo- 
nent which appears stronger than the core at this lower frequency. 
It could not be excluded, however, that the brightest component 
in Fig. ^ coincides with the core. This remains to be checked 
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with high-resolution VLBI imaging of J0839+0104 at frequencies 
higher than 1.6 GHz. 

Since the brightest component of the reference source, which 
is presumably not identical with the radio core, is in the phase cen- 
ter of our image, the relative position determination of the target 
source may be corrupted by a systematic error. Therefore we give 
a conservative error estimate of 8 mas for the position of SDSS 
0836+0054. This reflects the uncertainty in the identification of the 
reference point in the radio image of the phase-reference calibrator 
source J0839+0104. 

We searched for possible multiple images of SDSS 
0836+0054 that might have been caused by gravitational lensing. 
To this end, we investigated a large (4" x 4") field of view cen- 
tered around the a priori position of SDSS 0836+0054. To avoid 
distortion of the field of view due to frequency-average smearing, 
we averaged our data over only 8-MHz channels for imaging. We 
used the correlator output averaging time of 4 s which did not limit 
our field of view. No source of compact radio emission has been de- 
tected at the brightness level of ~ 100 ^ijy/beam (3 times the image 
noise), apart from what is seen in Fig. [2] Moreover, the total flux 
density in the CLEAN component model corresponding to the natu- 
rally weighted image of SDSS 0836+0054 (Fig.0 is 1.1 mjy. This 
is equal to the flux density in the VLA FIRST survey. Although 
the source total flux density recently measured with the VLA at 
1.4 GHz is 1.76 mjy (A. Petric et al., in preperation), the differ- 
ence may be attributed to source variability and, to a lesser extent, 
amplitude calibration uncertainties in our VLBI data as well. 

Our result suggests that most if not all of the radio emission 
of SDSS 0836+0054 is confined to a single compact object within 
an angular extent of ~ 10 mas. Thus we can conclude that the 
quasar is not multiply imaged by gravitational lensing at the level 
of brightness ratio ;£ 7. Outside of our field of view, there are no in- 
dications of multiple images neither in the VLA FIRST and NVSS 
surveys, nor in the Chandra X-ray image (e.g. lSchwartzl f2002). 
Most recentlv. lFan et aljiiooH) found that the high resolution (O'.'l) 
HST optical image of SDSS 0836+0054 is consistent with an un- 
resolved point source. However, it does not exclud e the pos sibility 
that its flux is magnified by a factor of up to ~ 2 IWvithe & Loebl 
2002b). Note that, based on the size of the ioinize d region around 
another SDSS quasar at z = 6.28, Haiman & Cen 1 2002) place an 
upper limit of 5 for the magnification factor in that case. 

Using the two-point radio spectral index of a = —0.9 (based 
on VLA total flux density measurements at 1 .4 and 5 GHz, A. Petric 
et al., in preperation; S cx u a , where S is the flux density and v 
is the frequency), the luminosity of SDSS 0836+0054 at 5 GHz 
(rest frame) is L 5 = 1.1 x 10 25 W Hz -1 sr" 1 . Togethe r with the 
estimated black hole mass of M hh = 4.8 x 1O 9 M iFan et ail 
1200 ll) . this value is consist ent with the L5 oc M^ h" relation found 
for radio quasars at z ^ 1 Ijarvis & McLurel200a and references 
therein). 



4 CONCLUSIONS 

SDSS 0836+0054, the most distant (z = 5.82) radio-emitting 
quasar known to date, is clearly detected with VLBI observations 
at 1.6 GHz using the EVN. A VLBI image with a dynamic range 
of 10:1 is made (Fig.|2j- The radio structure of the quasar at ~ 10- 
mas angular resolution is characterized by a compact component. 
Using phase-reference observations to a nearby calibrator, we de- 
rived the most accurate astrometric position of SDSS 0836+0054 
available at present. The accuracy of the coordinates is limited by 



the complex radio structure of the calibrator quasar PKS 0837+012, 
which we imaged for the first time at 1.6 GHz. The extremely dis- 
tant quasar SDSS 0836+0054 is unlikely to be multiply imaged by 
gravitational lensing. Magnification of its flux due to gravitational 
lensing or relativistic beaming cannot be ruled out. 
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